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In this paper we present an experimental approach that allows to deduce the important dynam- 
ical parameters of single sonoluminescing bubbles (pressure amplitude, ambient radius, radius-time 
curve) The technique is based on a few previously confirmed theoretical assumptions and requires 
the knowledge of quantities such as the amplitude of the electric excitation and the phase of the 
flashes in the acoustic period. These quantities are easily measurable by a digital oscilloscope, avoid- 
ing the cost of expensive lasers, or ultrafast cameras of previous methods. We show the technique 
on a particular example and compare the results with conventional Mie scattering. We find that 
within the experimental uncertainties these two techniques provide similar results. 
PACS numbers: 47.55.Dz, 43.25.+y, 78.60.Mq 



I. INTRODUCTION 

Single bubble sonoluminescence (SBSL) is a physical 
process where an oscillating gas bubble levitated acous- 
tically in a host liquid emits brief flashes of light in each 
period of the harmonic excitation [1]. In experiments the 
bubble levitation is done in resonators size of a jar, filled 
with liquid, and the acoustic excitation is accomplished 
by piezo-electric (PZT) transmitters [2]. The study of 
SBSL is motivated by the fact that through it one can 
study physical processes and the properties of matter at 
extreme conditions of high density, pressure and temper- 
ature. The variety of processes associated to the phe- 
nomenon include heat [3] and mass [4] transfer, chemical 
reactions [5] , shape oscillations [6] , chaos [7] , shock- wave 
[8] and light [9] emission to name a few. 

In order to validate theories of SBSL it is crucial to be 
able to measure several important parameters, such as 
the pressure amplitude near the bubble P^, the ambient 
radius i?o, or the radius-time curve R{t) in an acoustic 
period. Previous methods to measure these quantities are 
variations of Mie scattering technique [10], direct imag- 
ing of the bubble [11] and a method based on Doppler 
effect [12]. The common factor in these methods is that 
each of them requires a rather sophisticated and expen- 
sive experimental setup including lasers, precision optics, 
and high-speed cameras. Moreover each of these meth- 
ods are invasive^ i.e. they involve external action in the 
measuring process (e.g. laser or back-lighting), which 
may disturb the measurement of some other features of 
SBSL. The quantity Pa can also be measured directly by 
a needle hydrophone, however, because of the high cost 
and relatively low precision of this. Pa is usually deduced 
from fitting measured R{t) curves to the solutions of the 
Rayleigh-Plesset equation which describes the dynamics 
of the bubble's volumetric oscillations. 

These methods confirmed several theoretically pre- 
dicted features of SBSL, for example that stable light 
emitting bubbles obey diffusive and chemical stability 
[13-15]. In a previous study [16] we exploited this fact to 



develop a new technique to deduce P^, Ro or R{t). Here 
we elaborate on our method in more detail, present its 
advantages and limitations and compare its results with 
Mie scattering. 



II. DESCRIPTION OF THE FITTING 
TECHNIQUE 

Our tecfmique is based on the following assumptions: 
1. assumption Tfie dynamics of tlie bubble wall R{t) 
is well described by tfie Rayleigh-Plesset (RP) equation 
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where Pg is the uniform gas pressure inside the bubble, 
Pf = —Pa sm{ujt) is the forcing pressure with angular 
frequency cj, Pq is the ambient pressure valid during the 
measurements, and the remaining parameters are ma- 
terial constants of the host liquid, e.g. c is the speed of 
sound, p its density, and u is the kinematic viscosity. The 
gas pressure Pg can be related to R(t) through an equa- 
tion of state. We use a poly tropic van der Waals equation 
of state, modified to include the effects of surface tension 
(J and of vapor pressure P^, 
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Here a is the hard core van der Waals radius of the gas 
(for argon a = Ro/8.86), and 7 is the ratio of specific 
heats. In most of the acoustic period except the final 
stages of the collapse and after-bounces the gas can be 
considered isothermic [3], thus we use 7 = 1. 

2. assumption Bubbles emitting light in a stable fash- 
ion contain only inert gases, and are in stable diffusive 
equilibrium with the surrounding liquid [17]. The points 
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of diffusive equilibrium in the (Pa, -^o) space can be cal- 
culated from 



Ci/Co=^-^, {X) 



/o^" Rityxdt 
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where Pq* = 1 atm is the standard atmospheric pressure 
and Ta is the acoustic period. For air in water Ci is the 
concentration of argon set during the liquid preparation, 
Co is the tabulated dissolved air concentration of water 
under normal conditions. If the liquid preparation is done 
by degassing, then Ci can be calculated from Henry's law 
d = 0.0093Co- Pi /Pq*, where Pi is the partial pressure of 
air set during the preparation (air contains 0.93% of ar- 
gon). The diffusive equilibrium is stable where the curve 
in the (P^, Rq) space set by equation (3) is characterized 
by a positive slope. 3. assumption The acoustic pressure 
amplitude Pa at the bubble's position is directly propor- 
tional to the excitation voltage on the piezo transmitters 

Upzt- 



Pa = A-U, 



pzt 
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This assumption is reasonable, when the amplitude of 
displacement of the resonator walls can be considered 
small. For values typical in SBSL experiments P^ < 2 
bar this displacement is on the order of micrometers, 
thus equation (4) can be used safely. We define the di- 
mensionless phase of a flash ^ = tmm/Pa, as the time 
elapsed from the beginning of the acoustic period until 
the bubble reaches its minimum radius tmin normalized 
by the acoustic period. The phase ^ that can be calcu- 
lated numerically is the sum of the measured phase 
and of a constant parameter P, that accounts for possible 
phase-shift due to electronics. 
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The phase can be measured for instance by displaying the 
PMT signal of the flashes and the sinusoidal excitation 
voltage Upzt on an oscilloscope. 

The detailed procedure is the following. At an experi- 
mentally known dissolved gas concentration, liquid tem- 
perature T, ambient pressure and driving frequency one 
measures im{Upzt) at several excitation levels from the 
smallest possible light emission to the upper limit of SL 
[18]. After this equations (1), (2), (3) are solved system- 
atically in a wide range of the parameters Pa and Pq and 
the quantities ^, {Pg)^/PQ are extracted from the nu- 
meric R{t) data. From the resulting ^(Pa, Po, {^9)4^/^0) 
data one selects by interpolation a subset for which 
{Pg) ^ / P^ equals the experimentally known value of 
Ci/Co. This subset is the curve of diffusive equilibrium 
in the (Pa,Po5 space (see Fig. 1). The equilibrium is 
stable where dRo/dPa > or equivalently d^/dPa > 0. 
The final step is to fit the measured ^m{Upzt) data to 
the stable part of the diffusive equilibrium curve in the 
(Pa, ^) plane by adjusting the parameters A and B. This 



fit is highly constrained by the experimental fact that for 
a given dissolved gas concentration the lower limit of SL 
is linked to the smallest Pa on the stable diffusive equi- 
librium curve (see for instance Refs. [15,13,14]). After A 
and B is found the experimental data points can also be 
plotted in the (Pa, Po) plane and in principle R{t) curves 
and other dynamical parameters such as the expansion 
ratio Rmax/Ro are also determined. We stress here that 
only those parts of the numeric R{t) curves contain use- 
ful information about the size of the bubble, where the 
assumption of isothermicity 7=1 holds. These include 
the expansion phase, the maximum radius, and the ini- 
tial stages of the collapse, but exclude the region near 
the minimum radius and after bounces. 



III. EXPERIMENTAL APPARATUS AND THE 
DETAILS OF THE MEASUREMENT 

To test our method and to compare it with Mie scat- 
tering we performed the following experiment. In a de- 
gassing equipment [20] a distilled water sample was pre- 
pared with a dissolved air concentration 0.145Co at a 
temperature of T = 22^(7. This corresponds to an argon 
concentration of Ci = 0.00135Co. The water was trans- 
fered to the resonator [20] using gravity ffow. During 
the filling the gas pressure in the resonator was kept at 
the same level as in the degassing equipment, thus no air 
could diffuse in or out of the water. Such a care had to be 
taken because the precision of our method is guaranteed 
only if the gas concentration in the resonator is the same 
as set by the water preparation. After this the resonator 
was placed in a setup (Fig. 2) used for recording the SL 
flashes and the intensities of Mie scattered laser light. 

The digital oscilloscopes were triggered by the mon- 
itored driving signal Upzt- For each excitation level 
54 snapshots of Upzt^ scattered laser intensity and SL 
flashes were averaged, and the resulting traces were trans- 
fered to a PC. The measurement was done at an excita- 
tion frequency of 22700 Hz, the water temperature was 
T = 21.9 - 22.4^C7, and the external pressure Pq = 1017 
mbar. Figure 3(a) shows the fitting of the measured 
£,{Upzt) data to the calculated curve of stable diffusive 
equilibrium at Ci/Co = 0.00135. The error of ^ is deter- 
mined by the jitter in the phase of the flashes ~ 0.2 /is, 
thus ^ 0.00227 (the time resolution was 25 ns), while 
the error of the pressure amplitude APa = 0.01463 bar 
is calculated from 



APa = A • 2A5'. 
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The quantity A is the proportionality constant from 
equation (4), and AS = 0.001563F is the vertical res- 
olution of the digital scope corresponding to the 8-bit 
digitization. Both the zero-line and the maximum of the 
Upzt signal are known with a precision of AS, thus there 
is a factor of two in equation (6). Using the calculated 
^(Pa, Po) data and the values of ^ and Pa from the flt of 
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Fig. 3(a) one can plot the measured data points in the 
(Pa^Ro) plane (see the filled squares with error bars in 
Fig. 3(b)). The errors of ^ and Pa determine a range of 
Rq in the ^(Pa, ^o) data, which sets the error bar ARq. 



IV. COMPARISON WITH MIE SCATTERING 

The Mie scattering method is based on the detection 
of laser light scattered from the bubble, whose intensity 
under certain conditions is proportional to the square of 
the radius. More precisely the following relation holds 

R{tf = a{U{t)-Ubg), (7) 

where U{t) is the output signal of the PMT detecting the 
scattered intensity, and Ubg is the background scattered 
intensity in the absence of the bubble. The square-root 
of the recorded intensity is then fitted to a solution of the 
RP equation, and the desired parameters Pa, Ro, Rmax^ 
etc. are determined from the best fit. There are several 
difficulties of the technique that one must overcome. 

1. The scattered intensity has a very strong angle- 
dependence. This is reduced by choosing an appropriate 
angle and by use of a lens that averages out the light 
from different angles. In our measurement we used a 
setup similar to that of Barber et a/(1992) in [10]. 

2. The background scattered intensity is changing on 
a timescale of seconds because of small dust particles 
passing by near the position of the bubble due to a slow 
convection in the resonator. This can be handled partly 
by using an aperture and a laser beam as narrow as pos- 
sible, but can not be eliminated completely. 

3. Modification of the excitation level makes the bub- 
ble change its equilibrium position above the pressure 
anti-node as the averaged Bjerkness and buoyancy forces 
change [21]. The intensity of the laser light in the cross- 
section of the beam is a Gaussian, thus one always has 
to adjust the direction of the laser, in order to keep the 
bubble in the center of the beam. For this reason the pro- 
portionality constant a in equation (7) may be different 
for each excitation level depending on how precisely one 
can follow the bubble. Because of this difficulty we fitted 
the shape of the measured intensity signals for each exci- 
tation level, rather than using a proportionality constant 
found from a single fit. 

We used a numeric code that accomplished the fitting 
automatically in a user defined range of Pa and R^. The 
measured U (t) curve was transformed into a normalized 
positive signal u{t) = {U{t) - Uhg)/{U'' - Uhg), where t/* 
is the minimum of the negative U {t) which corresponds 
to the maximum radius of the bubble. Then a solution 
of the RP equation is generated at a given Pa and Ro, 
and the values of Rmax and tmin are extracted from the 
numeric radius- time curve. After this R{t) was trans- 
formed into r(t) = R{t - tmin + Kiin)l^max, where V^-^ 
is the time value of the end of the collapse in the mea- 
sured V(i) data. This way we got two normalized time 



series for which the phases corresponding to the mini- 
mum radius were identical. To be able to compare u(i) 
and r(t) these time series had to have the same number of 
elements, thus we interpolated r(t) at time values taken 
from u{t\ Finally the difference of the signals could be 
characterized by a single error parameter 

1=1 

where iV^^^ is the index of the minimum radius. The 
reason for calculating the error only until the end of the 
collapse is that in the after-bounce region the assumption 
of 7 = 1 used in equation (2) is not valid anymore. By 
this strategy for every measured U (t) curve we obtained 
the values of A{Pa,Ro) in a wide range of Pa and Rq 
and thus the dynamical parameters appropriate for the 
smallest error could be identified. Figure 4 shows the 
contour-plots of A for a particular case at two levels of 
the background. As can be seen the fits are best along a 
line in the {Pa, Rq) space, and the fitting is very sensitive 
to the value of U^g . 

In our measurements the precision of U^g was con- 
strained by the vertical resolution of the digital scope 
and the previously mentioned slow variance due to mov- 
ing dust particles to -0.0011 > Utg > -0.00225 (Volts). 
For comparison the values of C/* corresponding to the 
maximum radii were —0.05317.. — 0.07742 (Volts). In 
Fig. 4 the points within the contour A = 0.013 (which 
equals the noise level on the scattered intensity) deter- 
mine the errors of APa ~ 0.04 bar and ARq ^ 0.7 /im. If 
we also includes the uncertainty of Ubg then the final esti- 
mates are APa ~ 0.08 bar and ARq ~ 1 /im. These val- 
ues are consistent with the errors given by others [15,14]. 
The best fits with the bounding values of the background 
are shown in Fig. 3(b), where in the case of left and 
right triangles the background was Ubg = —0.0011, and 
Ubg = —0.00225 accordingly. The filled squares found 
from the fitting technique of sec. II in each case lie be- 
tween the best Mie fits confirming that these two meth- 
ods provide identical results within the experimental un- 
certainties. 

For further confirmation of the assumptions in sec. II 
we also fitted the Mie scattered data by a different strat- 
egy, where the background was allowed to vary between 
-0.0011 and -0.00225, and the R{t) data of the fits had 
to satisfy equation (3) exactly. The best fits of this strat- 
egy are shown as open circles in Fig. 3. As can be seen all 
the circles are within the error-bars found from the new 
technique. The quality of these Mie- fits are also shown 
individually in Figs. 5. 

For another independent test to confirm the assump- 
tions of the new technique we also analyzed measure- 
ments of Upzt and Mie scattered data of non-light- 
emitting "bouncing" bubbles. The open squares with 
error bars in Fig. 3(b) show the case where the pressure 
amplitudes were calculated from equation (4) using the 
value of A found from the fitting of the light emitting 
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bubbles. The ambient radii were determined afterward 
by best Mie fits at the prescribed Pa and by using the av- 
erage Ubg. Then the fitting of the Mie scattered data was 
also accomplished by using the same Ubg, but both the Pa 
and Rq parameters are allowed to vary, and the informa- 
tion of the Upzt data is not used. The best fits using this 
strategy are (shown as filled circles in Fig. 3(b)) in ex- 
cellent agreement with the best fits of the previous case. 
The phase diagram of Fig. 3(b) can also be compared to 
a measurement of reference [14], where the experimental 
conditions were quite close to ours (water was prepared 
with a dissolved air concentration of 0.14Co, and the ex- 
citation frequency was 20.6 KHz. The comparison with 
Fig. 4 of [14] shows a good quantitative agreement further 
justifying our method. 



V. CONCLUSION 

In this paper we presented an alternative method for 
the deduction of important dynamical parameters of 
SBSL, namely the pressure amplitude P^, the ambient ra- 
dius Rq, and the expansion ratio Rmax/Ro- The method 
is based on previously confirmed theoretical assumptions 
and the measurement of the phase of the flashes in the 
acoustic period together with the driving signal on the 
piezo transmitters, and the precise knowledge of the ex- 
perimental conditions, especially the dissolved gas con- 
centration in the liquid. We demonstrated the method 
by measuring the Pa and Rq parameters of light emit- 
ting bubbles for a given argon concentration. At the 
same time measurements of these parameters were also 
done by the conventional Mie scattering method, reveal- 
ing that within experimental uncertainties these two ap- 
proaches produced identical results. The phase diagram 
obtained is in good quantitative agreement with that of 
reference [14]. The limitation of the method is that unlike 
Mie scattering it can not be used to measure non-light- 
emitting bubbles, and also unable to measure unstable 
SL. However the important advantage of the method is 
its non-invasive chararcter, and thus the applicability to 
situations where the other methods are either imprac- 
tical or even potentially destructive to use as e.g. in 
single photon correlation measurements, or in measuring 
the parametric dependence of SBSL spectra. Moreover 
the technique does not require expensive apparatus other 
than a digital oscilloscope and a PMT, and the data anal- 
ysis is easily automatized. 
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FIG. 1. The curve of diffusive equilibrium in the (Pa,i?o,0 space, and its projections on the {Pa,Ro) and (Pa,0 planes for 
Ci/Co = 0.00135. Stable sonoluminescing bubbles follow the stable part of this curve (where dRo/dPa > 0, or d^/dPa > 0). 
The advantage of using the quantity ^ is that contrary to Rq it can be easily measured with high precision. 
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FIG. 2. The oscillating bubble was illuminated by a 30mW He-Ne laser, and the scattered laser light was focused by a 
lens through an aperture into the PMT (Hamamatsu R3478). Another PMT of the same kind detected the SL flashes. The 
signals of the PMTs were amplified and shaped by spectroscopy amplifiers (Ortec-model 451) and visualized together with the 
monitored electric signal of the piezo transmitters on digital oscilloscopes (HP-54616C 500MHz, HP-54600B lOOMHz) 
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FIG. 3. The fitting of the measured phases of the flashes to the numerically calculated curve (a), and the same data (filled 
squares with error bars) in the {Pa,Ro) plane (b). Also in (b) the left and right triangles correspond to best fits using Mie 
scattering and the bounding values for the background signal. The open circles stand for the best Mie fits assuming that 
equation (3) holds precisely, while the background was allowed to vary between the bounding values. The open squares indicate 
best Mie fits of non-light-emitting bubbles at Pa values found from the new technique, and finally the filled circles are best Mie 
fits where both Pa and Rq were fitted. The arrows indicate regions of shrinking and growing bubbles. 
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FIG. 5. The best fits of the normahzed u{t) and r{t)^ time series corresponding to the open circles in Fig. 3(b). The 
parameters of the fit are indicated in the corner of the figures. 
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